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Abstract 



The production of muons from heavy flavour decays is measured at forward rapidity in proton-proton 
collisions at y/s = 7 TeV collected with the ALICE experiment at the LHC The analysis is carried 
out on a data sample corresponding to an integrated luminosity L; nt = 16.5 nb _1 . The transverse 
momentum and rapidity differential production cross sections of muons from heavy flavour decays 
are measured in the rapidity range 2.5 < y < 4, over the transverse momentum range 2 < p t < 
12 GeV/c. The results are compared to predictions based on perturbative QCD calculations. 
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1 Introduction 



The study of heavy flavour (charm and beauty) production in proton-proton collisions at LHC (Large 
Hadron Collider) energies provides an important test of perturbative QCD (pQCD) calculations 0] in 
a new energy domain, where unprecedented small Bjorken-x (momentum fraction) values are probed. In 
the rapidity region 2.5 < y < 4, charm (beauty) production at yfs = 7 TeV is expected to be sensitive to 
x values down to about 6 • 10~ 6 (2 • 10~ 5 ). Important progress has been achieved in the understanding of 
heavy flavour production at lower energies. In earlier measurements, the beauty production cross section 
in pp collisions at \fs =1.8 TeV measured by the CDF and DO experiments [3,4] at the FNAL Tevatron, 
was found to be higher than Next to Leading Order (NLO) pQCD predictions HI . More recent results 
from the CDF Collaboration Q, for pp collisions at y/s = 1.96 TeV, are described well by Fixed Order 
Next-to-Leading Log (FONLL) OH and NLO [8] pQCD calculations. The charm production cross 
section measured at the FNAL Tevatron [9] is also well reproduced by FONLL [10] and GM-VFN iTTTt 
calculations within experimental and theoretical uncertainties, although at the upper limit of the calcu- 
lations. The PHENIX and STAR Collaborations HH [131 at the RHIC (Relativistic Heavy Ion Collider) 
measured the production of muons and electrons from heavy flavour decays in pp collisions at y/s = 0.2 
TeV The upper limit of FONLL pQCD calculations 041 is consistent with the measurement of electrons 
from heavy flavour decays in the mid-rapidity region, while in the forward rapidity region the produc- 
tion of muons from heavy flavour decays is underestimated by the model calculations. Furthermore, at 
LHC energies, the ATLAS 03), LHCb 02) and CMS 07JQ2) Collaborations reported on the measure- 
ment of beauty production in pp collisions at y/s = 7 TeV. The results are consistent with NLO pQCD 
calculations within uncertainties. A similar agreement with FONLL calculations is also observed for 
mid-rapidity electrons and muons from heavy flavour decays, measured by the ATLAS experiment lfl9l 
in pp collisions at yfs = 7 TeV. In this respect, it is particularly interesting to perform the measurement 
of heavy flavour decay muon production in the forward rapidity region at the LHC and compare it with 
theoretical models. 

The investigation of heavy flavour production in pp collisions also constitutes an essential baseline for the 
corresponding measurements in heavy-ion collisions. In the latter, heavy quarks are produced at early 
stages of the collision and then experience the full evolution of the extremely hot and dense, strongly 
interacting medium ll20ll2TI . The modification of the heavy flavour transverse momentum distributions 
measured in heavy ion collisions with respect to those measured in pp collisions is considered as a 
sensitive probe of this medium |[22ll23l . 

Finally, the study of heavy flavour production is also important for the understanding of quarkonium 
production, both in pp, p-nucleus and nucleus-nucleus collisions If20l [2T1 . 

The ALICE experiment ll24l measures the heavy flavour production at mid-rapidity through the semi- 
electronic decay channel G51 and in a more direct way through the hadronic D-meson decay chan- 
nel J26), and at forward rapidity through the semi-muonic decay channel. In this paper, we present the 
measurement of differential production cross sections of muons from heavy flavour decays in the rapidity 
range 2.5 < y < 4 and transverse momentum range 2 < p t < 12 GeV/c, with the ALICE muon spectrom- 
eter ll24l . in pp collisions at y/s = 7 TeV. The results are compared to FONLL pQCD calculations |[2ll27). 

The paper is organized as follows. Section 2 consists of an overview of the ALICE experiment with an 
emphasis on the muon spectrometer and a description of data taking conditions. Section 3 is devoted 
to the analysis strategy: event and track selection, background subtraction, corrections, normalization 
and determination of systematic uncertainties. Section 4 addresses the experimental results: p t - and 
^-differential production cross sections of muons from heavy flavour decays at forward rapidity, and 
comparisons to FONLL pQCD predictions. Conclusions are given in Section 5. 
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2 The ALICE experiment and data taking conditions 

A detailed description of the ALICE detector can be found in ll24l . The apparatus consists of two 
main parts: a central barrel (pseudo-rapidity coverage: \r} \ < 0.9) placed in a large solenoidal mag- 
net (B = 0.5 T), which measures hadrons, electrons and photons, and a muon spectrometer (—4 < T] < 
— 2.5^]). Several smaller detectors for global event characterization and triggering are located in the for- 
ward and backward pseudo-rapidity regions. Amongst those, the VZERO detector is used for triggering 
purposes and in the offline rejection of beam-induced background events. It is composed of two scintilla- 
tor arrays placed at each side of the interaction point and covering 2.8 < v\ < 5.1 and —3.7 < r\ < — 1.7. 
The central barrel detector used in this work for the interaction vertex measurement is the Silicon Pixel 
Detector (SPD), the innermost part of the Inner Tracking System (ITS). The SPD consists of two cylin- 
drical layers of silicon pixels covering \r\ \ < 2.0 and \r\\ < 1 .4 for the inner and outer layer, respectively. 
The SPD is also used in the trigger logic. 

The muon spectrometer detects muons with momentum larger than 4 GeV/c and is composed of two 
absorbers, a dipole magnet providing a field integral of 3 Tm, and tracking and trigger chambers. A 
passive front absorber of 10 interaction lengths (Ai), made of carbon, concrete and steel, is designed to 
reduce the contribution of hadrons, photons, electrons and muons from light hadron decays. A small 
angle beam shield (6 < 2°), made of tungsten, lead and steel, protects the muon spectrometer against 
secondary particles produced by the interaction of large-T] primary particles in the beam pipe. Tracking 
is performed by means of five tracking stations, each composed of two planes of Cathode Pad Chambers. 
Stations 1 and 2 (4 and 5) are located upstream (downstream) of the dipole magnet, while station 3 is 
embedded inside the dipole magnet. The intrinsic spatial resolution of the tracking chambers is better 
than 100 /im. Two stations of trigger chambers equipped with two planes of Resistive Plate Chambers 
each are located downstream of the tracking system, behind a 1.2 m thick iron wall of 7.2 At. The latter 
absorbs most of the hadrons that punch through the front absorber, secondary hadrons produced inside 
the front absorber and escaping it and low momentum muons (p < 4 GeV/c). The spatial resolution of 
the trigger chambers is better than 1 cm and the time resolution is about 2 ns. Details concerning track 
reconstruction can be found in |[28l l29l . 

The results presented in this publication are based on the analysis of a sample of pp collisions at 
s/s = 7 TeV collected in 2010, corresponding to an integrated luminosity of 16.5 nb _1 . 

The data sample consists of minimum bias trigger events (MB) and muon trigger events (/X-MB), the 
latter requiring, in addition to the MB trigger conditions, the presence of one muon above a transverse 
momentum (p t ) threshold that reaches the muon trigger system. The MB trigger is defined as a logical 
OR between the requirement of at least one hit in the SPD and a hit in one of the two VZERO scintillator 
arrays. It also asks for a coincidence between the signals from the two beam counters, one on each side 
of the interaction point, indicating the passage of bunches. This corresponds to at least one charged 
particle in 8 units of pseudo-rapidity. The logic of the jU-MB trigger requires hits in at least three (out of 
four possible) trigger chamber planes. The estimate of the muon transverse momentum is based on the 
deviation of the measured track with respect to a straight line coming from the interaction point, in the 
bending plane (plane measuring the position along the direction perpendicular to the magnetic field). By 
applying a cut on this deviation, tracks above a given p t threshold are selected. The p t threshold allows 
the rejection of soft background muons mainly coming from pion and kaon decays, and also to limit 
the muon trigger rate when high luminosities are delivered at the interaction point. In the considered 
data taking period, the p t trigger threshold was set to its minimum value of about 0.5 GeV/c and the 
corresponding muon trigger rate varied between about 40 and 150 Hz. The instantaneous luminosity 
at the ALICE interaction point was limited to 0.6-1.2- 10 29 cm s by displacing the beams in the 

'The muon spectrometer covers a negative pseudo-rapidity range in the ALICE reference frame, r/ and y variables are 
identical for muons in the acceptance of the muon spectrometer, and in pp collisions the physics results are symmetric with 
respect to 7] (y) = 0. They will be presented as a function of y, with positive values. 
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transverse plane by 3.8 times the r.m.s of their transverse profile. In this way, the probability to have 
multiple MB interactions in the same bunch crossing is kept below 2.5%. 

The alignment of the tracking chambers, a crucial step for the single muon analysis, was carried out 
using the MILLEPEDE package IT301 . by analyzing tracks without magnetic field in the dipole and 
solenoidal magnet. The corresponding resolution is about 300 /xm in the bending plane, for tracks with 
p t > 2 GeV/c. With such alignment precision, the relative momentum resolution of reconstructed tracks 
ranges between about 1% at a momentum of 20 GeV/c and 4% at 100 GeV/c. 

3 Data analysis 

The single muon analysis was carried out with muon trigger events while, as will be discussed in Sec- 
tion 3.4, minimum bias trigger events were used to convert differential muon yields into differential cross 
sections. The identification of muons from charm and beauty decays in the forward region is based on the 
p t distribution of reconstructed tracks. Three main background contributions must be subtracted and/or 
rejected: 

- decay muons: muons from the decay of primary light hadrons including pions and kaons (the main 
contribution) and other meson and baryon decays (such as J/y and low mass resonances TJ, p, CO 
and </>); 

- secondary muons: muons from secondary light hadron decays produced inside the front absorber; 

- punch-through hadrons and secondary hadrons escaping the front absorber and crossing the track- 
ing chambers, which are wrongly reconstructed as muons. 

A Monte Carlo simulation based on the GEANT3 transport code EQISll and using the PYTHIA 6.4.21 
event generator 1331 l34l (tune Perugia-0 11351 ) was performed to obtain the p t distributions of these 
different contributions. They are displayed in Fig. [T] after all the selection cuts discussed in Section 
3.1 were applied. After cuts, the component of muons from heavy flavour decays prevails over the 
background contribution for p t > 4 GeV/c. The simulation results indicate that the hadronic background 
and the contribution of fake tracks (tracks which are not associated to one single particle crossing the 
whole spectrometer) are negligible. The component of muons from and Z° decays, which dominates 
in the p t range 30-40 GeV/c 11361 [T9l . is not considered in this analysis. This contribution is negligible 
in the p t range of interest 2-12 GeV/c. 

3.1 Data sample: event and track selection 

The data sample used in the physics analysis amounts to 1.3T0 7 jlt-MB trigger events. These selected 
events satisfied the quality criteria on detector conditions during data taking and the analysis quality cri- 
teria, which reduced the beam-induced background. This was achieved by using the timing information 
from the VZERO and by exploiting the correlation between the number of hits and track segments in 
the SPD. The accepted events have at least one interaction vertex reconstructed from hits correlation in 
the two SPD layers. The corresponding total number of tracks reconstructed in the muon spectrometer 
is 7.8-10 . Various selection cuts were applied in order to reduce the background contributions in the 
data sample. Tracks were required to be reconstructed in the geometrical acceptance of the muon spec- 
trometer, with — 4 < r\ < —2.5 and 171° < a t, s < 178°, 6 a \, s being the track polar angle measured at 
the end of the absorber. These two cuts reject about 9% of tracks. Then, the track candidate measured 
in the muon tracking chambers was required to be matched with the corresponding one measured in the 
trigger chambers. This results in a very effective rejection of the hadronic component that is absorbed 
in the iron wall. This condition is fulfilled for a large fraction of reconstructed tracks since the anal- 
ysis concerns jii-MB trigger events. The fraction of reconstructed tracks that are not matched with a 
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Fig. 1: Transverse momentum distribution of reconstructed tracks in the muon spectrometer after all selection 
cuts were applied (see section 3.1 for details). The distributions were obtained from a PYTHIA 11331 l34l (tune 
Perugia-0 [35 1) simulation of pp collisions at y/s = 7 TeV. The main sources are indicated in the figure. 



corresponding one in the trigger system is about 5%. For comparison, in MB collisions this fraction is 
about 64%. Furthermore, the correlation between momentum and Distance of Closest Approach (DCA, 
distance between the extrapolated muon track and the interaction vertex, in the plane perpendicular to 
the beam direction and containing the vertex) was used to remove remaining beam-induced background 
tracks which do not point to the interaction vertex. Indeed, due to the multiple scattering in the front 
absorber, the DCA distribution of tracks coming from the interaction vertex is expected to be described 
by a Gaussian function whose width depends on the absorber material and is proportional to I /p. The 
beam-induced background does not follow this trend and can be rejected by applying a cut on p x DCA 
at 5a, where a is extracted from a Gaussian fit to the px DCA distribution measured in two regions in 
a bs, corresponding to different materials in the front absorber. This cut removes 0.4% of tracks, mainly 
located in the high p t range (in the region p t > 4 GeV/c, this condition rejects about 13% of tracks). 
After these cuts, the data sample consists of 6.67-10 6 muon candidates. 

The measurement of the heavy flavour decay muon production is performed in the region p t > 2 GeV/c 
where the contribution of secondary muons is expected to be small (about 3% of the total muon yield, 
see Fig. [T}. In such a p t region the main background component consists of decay muons and amounts 
to about 25% of the total yield (see Fig.[T]). 

3.2 Subtraction of the background contribution of decay muons 

The subtraction of the background component from decay muons (muons from primary pion and kaon 
decays, mainly) is based on simulations, using PYTHIA 6.4.21 Il33ll34l (tune Perugia-0 ll35l ) and PHO- 
JET 1.12 |37l as event generators. In order to avoid fluctuations due to the lack of statistics in the high 
p t region in the Monte Carlo generators, the reconstructed p t distribution of decay muons, obtained after 
all selection cuts are applied (Section 3.1), is fitted using 

di t = \ 7 pJ+W , (1) 

where a, b and c are free parameters. The fits are performed in five rapidity intervals, in the region 
2.5 < y < 4. The normalization is done assuming that the fraction of decay muons in the data is the 
same as the one in the simulations, in the region where this component is dominant (p t < 1 GeV/c). 
Finally, the (fitted) p t distribution is subtracted from the measured muon p t distribution. The subtracted 
p t distribution is the mean of the p t distributions from the PYTHIA and PHOJET event generators. 
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Table 1: Systematic uncertainties introduced by the procedure used for the subtraction of decay muons. MC and 
transport refer to the systematic uncertainty due to model input and transport code, respectively. See the text for 
details. 
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The total systematic uncertainty due to this procedure includes contributions from the model input and 
the transport code (GEANT3 0T1 1321 "). The former takes into account the shape and normalization of 
the p t distribution of decay muons, and the observed difference in the K 1 * 1 / ratio as a function of p t in 
the mid-rapidity region [38] between ALICE data and simulations. The results show that both PYTHIA 
(tune Perugia-0) and PHOJET underestimate this ratio by about 20%. The corresponding uncertainty due 
to this difference between data and simulations is propagated to the muon yield in the forward rapidity 
region. The effect of the transport code is estimated by varying the yield of secondary muons within 
100% in such a way to provide a conservative estimate of the systematic uncertainty on the secondary 
particle production in the front absorber. The systematic uncertainty from the model input varies from 
about 7% to 2% as y increases from 2.5 to 4, independently of p t , while the one from the transport code 
depends both on y and p t and ranges from 4% (3.7 < y < 4) to a maximum of 34% (p t = 2 GeV/c and 
2.5 < y < 2.8). The corresponding values of these systematic uncertainties as a function of p t and y are 
summarized in Table [T] They are added in quadrature in the following. 

3.3 Corrections 

The extracted yields of muons from heavy flavour decays are corrected for acceptance, reconstruction 
and trigger efficiencies by means of a simulation modelling the response of the muon spectrometer. 
The procedure is based on the generation of a large sample of muons from beauty decays by using a 
parameterization of NLO pQCD calculations ||29l . The tracking efficiency takes into account the status 
of each electronic channel and the residual mis-alignment of detection elements. The evolution of the 
tracking efficiency over time is controlled by weighting the response of electronic channels as a function 
of time. The typical value of muon tracking efficiency is about 93%. The efficiencies of the muon trigger 
chambers are obtained directly from data [28 ] and employed in the simulations. The typical value of such 
efficiencies is about 96%. Figure [2] shows the resulting acceptance and efficiency (A x e) as a function 
of generated p t . The global Axe increases significantly up to about 2 GeV/c and tends to saturate at a 
value close to 90%. 

The systematic uncertainty corresponding to the sensitivity of A x e on the input p t and y distributions 
was estimated by comparing the results with those from a simulation using muons from charm decays. 
This amounted to less than 1% and was neglected. The accuracy in the detector modelling introduces a 
systematic uncertainty estimated to be 5%, by comparing the values of trigger and tracking efficiencies 
extracted from data and simulations [28]. 

The distortion of the measured p t distribution, dominated in the high p t region by the effect of residual 
mis-alignment, is also corrected for by introducing in the simulation a residual mis-alignment of the 
same order of magnitude as in the data. However, this residual mis-alignment is generated randomly. 
A p t dependent relative systematic uncertainty on the muon yield of 1% x p t (in GeV/c) is considered 
in order to take into account the differences between the real (unknown) residual mis-alignment and the 
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simulated one. This is a conservative value determined by comparing the reconstructed p t distribution 
with or without including the residual mis-alignment. 
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Fig. 2: Acceptance x efficiency as a function of generated p t , obtained from a simulation of muons from beauty 
decays. 



3.4 Production cross section normalization 

The differential production cross section is obtained by normalizing the corrected yields of muons from 
heavy flavour decays to the integrated luminosity. Since the yields have been extracted using /i-MB 
trigger events, the differential production cross section is calculated according to: 

d 2 <V^HF _ d 2 ^±^HF A^MB a MB 

dp t dy d Pt dy N^_ MB n mb 

where: 

- — dpidy" 17 * s t ' ie P { an( * y differential yield of muons from heavy flavour decays; 

„± „± 

- A^ B and N^_ MB are the numbers of reconstructed tracks that satisfy the analysis cuts in MB and 
jU-MB trigger events, respectively; 

- jVmb is the number of minimum bias collisions corrected as a function of time by the probability to 
have multiple MB interactions in a single bunch crossing, and Omb is the corresponding measured 
minimum bias cross section. 



Cmb is derived from the C7vzero-and cross section [39] measured with the van der Meer scan method |40l . 
The VZERO — AND condition is defined as a logical AND between signals in the two VZERO scintil- 
lator arrays. Such a combination allows one to reduce the sensitivity to beam-induced background. The 
<7vzero-and/0mb ratio is the fraction of minimum bias events where the VZERO — AND condition is 
fulfilled. Its value is 0.87 and it remains stable within 1% over the analyzed data sample. This gives Omb 
= 62.5 ± 2.2 (syst.) mb. The statistical uncertainty is negligible, while the 3.5% systematic uncertainty 
is mainly due to the uncertainty on the beam intensities BTI and on the analysis procedure related to the 
van der Meer scan of the VZERO — AND signal. Other effects, such as oscillation in the ratio between 
MB and VZERO - AND counts, contribute less than 1%. 

3.5 Summary of systematic uncertainties 

The systematic uncertainty on the measurements of the p t and y differential production cross sections 
of muons from heavy flavour decays accounts for the following contributions discussed in the previous 
sections: 
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- background subtraction: from about 5% (3.7 < y < 4) to a maximum of 35% (2.5 < y < 2.8, 
Pt = 2 GeV/c), see Section 3.2 and Table[TJ 

- detector response: 5% (Section 3.3); 

- residual mis-alignment: 1% x p t (Section 3.3); 

- luminosity measurement: 3.5% (Section 3.4). 

The resulting systematic uncertainty, in the rapidity region 2.5 < y < 4, varies between 8-14% (the 3.5% 
systematic uncertainty on the normalization is not included). 

4 Results and model comparisons 



The measured differential production cross sections of muons from heavy flavour decays as a function 
of p t in the rapidity region 2.5 < y < 4 and as a function of y in the range 2 < p t < 12 GeV/c are 
displayed in Fig. [3] (circles), left and right panels, respectively. The error bars (which are smaller than 




p (GeV/c) ... ^ 

Fig. 3: Left: /^-differential production cross section of muons from heavy flavour decays in the rapidity range 
2.5 < y < 4. Right: y-differential production cross section of muons from heavy flavour decays, in the range 
2 < pt < 12 GeV/c. In both panels, the error bars (empty boxes) represent the statistical (systematic) uncertainties. 
A 3.5% normalization uncertainty is not shown. The solid curves are FONLL calculations and the bands display 
the theoretical systematic uncertainties. Also shown, are the FONLL calculations and systematic theoretical un- 
certainties for muons from charm (long dashed curves) and beauty (dashed curves) decays. The lower panels show 
the corresponding ratios between data and FONLL calculations. 

symbols in most of the p t and y bins) represent the statistical uncertainties. The boxes correspond to the 
systematic uncertainties. The systematic uncertainty on Omb is not included in the boxes. The results are 
compared to FONLL predictions El |2Vj (black curve and shaded band for the systematic uncertainty). 
The central values of FONLL calculations use CTEQ6.6 [42] parton distribution functions, a charm 
quark mass (ra c ) of 1.5 GeV/c 2 , a beauty quark mass (rat,) of 4.75 GeV/c 2 and the renormalization 

(jUr) and factorization (jUf) QCD scales such that jUr/jUo = Mf/Mo = 1 (Mo = ^t,q = J Pt+ m \)- The 
theoretical uncertainties correspond to the variation of charm and beauty quark masses in the ranges 
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1.3 < m c < 1.7 GeV/c 2 and 4.5 <m b < 5.0 GeV/c 2 , and QCD scales in the ranges 0.5 < jUr/jUq < 2 and 
0.5 < jUf/jUo < 2 with the constraint 0.5 < jIIf/Mr < 2. The FONLL predictions for muons from beauty 
decays include the components of muons coming from direct b-hadron decays and from b-hadron decays 
via c-hadron decays (e.g. B — >• D — > }X channel). The uncertainty band is the envelope of the resulting 
cross sections. The ratios between data and FONLL predictions are shown in the bottom panels. A good 
description of the data is observed within uncertainties, for both the p t distribution (up to 12 GeV/c) 
and the y distribution (in the p t range from 2 to 12 GeV/c). The measured production cross sections 
are systematically larger than the central values of the model predictions. The ratio of data over central 
values of FONLL calculations as a function of p t and y is about 1.3 over the whole p t and y ranges. 
This is consistent with the ALICE measurements of the /^-differential production cross sections of D 
mesons ll26l in the central rapidity region. The CMS and ATLAS Collaborations made complementary 
measurements of the heavy flavour production, with electrons and/or muons measured at mid-rapidity in 
pp collisions at y/s = 7 TeV H8lH9ll . The production of muons from beauty decays, measured by the 
CMS Collaboration in |tj| < 2.1 and at high p t (p t > 6 GeV/c), exhibits a similar agreement with NLO 
pQCD calculations within uncertainties: the data points lie in the upper limit of the model predictions. 
The results from the ATLAS Collaboration concerning the production of muons and electrons from heavy 
flavour decays in |t]| < 2.0 (excluding 1.37 < |t]| < 1.52) and in the region 7 < p t < 27 GeV/c are also 
consistent with FONLL calculations. 

The theoretical charm and beauty components are also displayed in Fig. [3] According to these predic- 
tions, the muon contribution from beauty decays is expected to dominate in the range p t > 6 GeV/c. In 
this region, it represents about 62% of the heavy flavour decay muon cross section. 

A similar comparison between data and FONLL calculations was performed in five rapidity intervals 
from y = 2.5 to y = 4 (Fig. |4| upper panels). The corresponding ratio of data over FONLL predictions is 
depicted in the lower panels of Fig.|4] The model calculations provide an overall good description of the 
data up to p t = 12 GeV/c in all rapidity intervals, within experimental and theoretical uncertainties. 

5 Conclusions 

We have presented measurements of the differential production cross sections of muons from heavy 
flavour decays in the rapidity range 2.5 < y < 4 and transverse momentum range 2 < p t < 12 GeV/c, in 
pp collisions at y/s = 7 TeV with the ALICE experiment. The FONLL pQCD calculations are in good 
agreement with data within experimental and theoretical uncertainties, although the data are close to the 
upper limit of the model calculations. Both the p t and y dependence of the heavy flavour decay muon 
production cross section is well described by the model predictions. The results provide an important 
baseline for the study of heavy quark medium effects in nucleus-nucleus collisions. 
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Fig. 4: Upper panel: /^-differential production cross section of muons from heavy flavour decays in five rapidity 
regions mentioned in the figures. The error bars (empty boxes) represent the statistical (systematic) uncertainties. 
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the theoretical systematic uncertainty. Lower panel: ratio between data and FONLL calculations. 
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